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Abstract — The eigenvalue spectrum of the adjacency matrix 
of a network is closely related to the behavior of many dynami- 
cal processes run over the network. In the field of robotics, this 
spectrum has important implications in many problems that 
require some form of distributed coordination within a team of 
robots. In this paper, we propose a continuous-time control 
scheme that modifies the structure of a position-dependent 
network of mobile robots so that it achieves a desired set of 
adjacency eigenvalues. For this, we employ a novel abstraction 
of the eigenvalue spectrum by means of the adjacency matrix 
spectral moments. Since the eigenvalue spectrum is uniquely 
determined by its spectral moments, this abstraction provides 
a way to indirectly control the eigenvalues of the network. Our 
construction is based on artificial potentials that capture the 
distance of the network's spectral moments to their desired 
values. Minimization of these potentials is via a gradient descent 
closed-loop system that, under certain convexity assumptions, 
ensures convergence of the network topology to one with the 
desired set of moments and, therefore, eigenvalues. We illustrate 
our approach in nontrivial computer simulations. 

I. Introduction 

A wide variety of coordinated tasks performed by teams 
of mobile robots critically rely on the topology of the un- 
derlying communication network and its spectral properties. 
Examples include, cooperative manipulation [l]-[3], surveil- 
lance and coverage [4]-[6], distributed averaging [7]-[9], for- 
mation control [ 10]— [ 12], flocking [13, 14], and multi-robot 
placement [15]— [17], that all require some form of network 
connectivity, structure and, oftentimes, spectral properties. In 
this paper we address the problem of controlling a network 
of mobile robots to a topology with a desired eigenvalue 
spectrum. This effort is a first step towards the design of 
controllers that allow robots to perform their assigned tasks, 
while optimizing coordination within the team. 

The eigenvalue spectra of a network provide valuable 
information regarding the behavior of many dynamical pro- 
cesses running within the network [18]. For example, the 
eigenvalue spectra of the Laplacian and adjacency matrices 
of a graph affects the mixing speed of Markov chains [19], 
the stability of synchronization of a network of nonlinear 
oscillators [20, 21], the spreading of a virus in a network 
[22, 23], as well as the dynamical behavior of many decen- 
tralized network algorithms [24]. Similarly, the second small- 
est eigenvalue of the Laplacian matrix (also called spectral 
gap) is broadly considered a critical parameter that influences 
the stability and robustness properties of dynamical systems 
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that are implemented over information networks [25, 26]. 
Optimization of the spectral gap has been studied both in a 
centralized [27]-[29] and decentralized context [30]. 

In this paper, we propose a novel framework to control 
the structure of a network of mobile robots to achieve a 
desired eigenvalue spectrum. In particular, we focus on the 
spectrum of the weighted adjacency matrix of the network, 
with weights that are (decreasing) functions of the inter- 
robot distances. This construction is relevant, for example, 
in modeling the signal strength in wireless communication 
networks. Although our framework performs well with dif- 
ferent distance metrics, in this paper, we focus on the t\ 
norm (Manhattan distance) primarily for analytical reasons, 
since its composition with convex functions preserves their 
convexity. Additionally, this metric has potential applications 
in indoor navigation where the presence of obstacles forces 
the signals to propagate in grid-like environments. 

We employ a novel abstraction of the eigenvalue spectrum 
in terms of the associated spectral moments, and define 
artificial potentials that capture the distance between the 
network's spectral moments and their desired values. These 
potentials are minimized via a gradient descent algorithm, 
for which we show convergence to the globally optimal mo- 
ments. Since the eigenvalue spectrum is uniquely determined 
by the associated spectral moments, our approach provides 
a way to indirectly control a network's eigenvalues. This 
work is related to [31], which addresses a similar problem 
for static robots in discrete environments. This formulation, 
however, is more appropriate for robotics applications, where 
communication depends continuously on the robot motion. 

The rest of this paper is organized as follows. In Section II- 
|A"1 we introduce some graph-theoretical notation and useful 
results. In Section |II-B| we formulate the control problem 
under consideration. We introduce an artificial potential and 
derive the associated motion controllers in Section IIII-AI 



and discuss convergence of our approach in Section III-B 



Finally, in Section IV we illustrate our approach with several 
computer simulations. 

II. Preliminaries & Problem Definition 

A. Notation and Preliminaries 

In this section we introduce some nomenclature and results 
needed in our exposition. Let ^ = {~f,S,W) denote a 
weighted undirected graph, with ~f = [n] being a set of n 
nodes, S C ~f x "V a set of e undirected edges, and f ei f + 
a set of weights associated to the edges. If {i,j} £ S we call 
nodes i and j adjacent (or neighbors), which we denote by 
i ~ j. In this paper, we consider graphs without self-loops, 
i.e., {z, i] £ § for all i £ V . We denote by aij = aji £ R+ the 



weight associated with edge {i, j} £ §, and assume that au = 
for {i,j} ^ £ '. We define a walk w of length k from node 
v'o to Vk to be an ordered sequence of nodes (vo,Vi, vj) 
such that Vj ~ v,+i for / = 0, 1, ...,k— 1. We define the weight 
a w of a walk w = (v , vi, v*) as a w = nf=o Ovj,v /+ i • 

A weighted graph can be algebraically represented via its 
weighted adjacency matrix, defined as the nx n symmetric 
matrix Ay — [ay], where ay is the weight of edge {i,j}. 
In this paper we are particularly interested in the spectral 
properties of Ay. Since Ay is a real symmetric matrix, it has 
a set of real eigenvalues Ai > X2 > ... > A„. The powers of 
the adjacency matrix Ay can be related to walks in C S: 

Lemma 1: Let be a weighted undirected graph with no 
self-loops. The (i,j)-th entry of A|, can be written in terms 
of walks in as follows: 



WGW; 



where W/,- is the set of all closed walks of length k from 
node i to j in the complete graph K,^ 

Proof: The above result comes directly from an alge- 
braic expansion of [A^] y . First, notice that 



A k 
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Using the above rule in a simple recursion, we can expand 
the entries of the decreasing powers of Ay, to obtain 
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which is the statement of our lemma. ■ 
The following result will be useful in our derivations: 
Lemma 2: Let be a weighted undirected graph with no 

self-loops. Then, we have that 

d\x{A%) 



da; 



= 2k 



(1) 



Proof: First, notice that for any two matrices S = [sy] 
and B = [bij], we have that tr(SB) = Y*ij s ijbji- Consider S 
and B to be symmetric matrices. Then, we can write S = U + 
L, where U and L are upper and lower triangular matrices, 
respectively, with L = U T . Let i < j (the same holds for 
j < i), hence 
d 



dst, 



tr(£B) = 
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'The complete graph ST„, is the undirected graph with n nodes in which 
every pair of distinct vertices is connected by a unique edge. 



Then, we have that 
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which completes the proof. ■ 

B. Problem Definition 

Consider a group of n mobile robots and define by 
Xi(t) £ M. d the position of robot i at time t > 0. Let x — 
[x\ . . . x%\ € M. dn denote the stacked column vector of 
all robot positions, so that x,> is the r-th coordinate of the z'-th 
robot position. We assume that we can control the position 
of the robots by the simple kinematic law 



Diln 



(2) 



where m,- (f ) £ W 1 is the control input applied to robot i, and 
Uj r (t) is the r-th coordinate of Hj(f). 

For a given set of robot positions, {jc,}" =1 , we define 
the weighted adjacency matrix of the network of robots as 

A(x) — [(tij(x)] with 



aij(x) 



(3) 



where c > is a constant and z £ {1,2} denotes the l\ or £2 
norm. Notice that A(x) is a symmetric matrix with position- 
dependent real eigenvalues, {A, (x)}" =1 . We define, further, 
the fc-th spectral moment of A (x) by 



m k (x)±-£ j Xf(x)= 1 -tr\A k (x) 
",ti n L 



(4) 



for 1 < k < n, where the last equality follows from diag- 
onalizing A (x). For any finite network with n nodes, its 
eigenvalue spectrum is uniquely defined by the sequence of 
n moments (mi,m2, ...,m„). Therefore, we can simultane- 
ously control the whole set of eigenvalues of a network by 
controlling the first n spectral moments of A(x). This gives 
rise to the following problem that we aim to address: 

Problem 3 (Control of spectral moments): Let { m t}1 = \ 
denote a desired set of spectral moments. Design control 
laws m; for all robots i = 1,. . . ,n so that the adjacency matrix 
A(x) of the position-dependent robot network has spectral 
moments that satisfy m^x) — > m\ for all k — l,...,n. 

Controlling the eigenvalue spectrum of the adjacency 
matrix of a network is particularly important, since it is 
related to the behavior of interesting network dynamical 
properties [18]. In Section [Hi] we propose a gradient descent 
algorithm to address Problem [3] and discuss its convergence 
properties. In Section |IV] we illustrate our approach in 
numerical simulations. 



III. Control of Spectral Moments 

A. Controller Design 

Assume a given sequence of desired spectral moments 
{«tt}5_ j and define the cost function 

1 



/(*)- E ^.( m k( x )- m k) 



k=l 



(5) 



where m k (x) is defined in (j4j). Define, further, the gradient 
descent control law 



-V*/(x) 



(6) 



Then, we can show the following result: 

Lemma 4: Let the adjacency matrix A=A(x) be defined 
as in ^ and assume that z = 1 (ti norm). Then, an explicit 
expression for the entries of u is given by 



k=l 



trA K 



■m. 



(AoS,-)A 
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where S r (x) — [sjj] with sy =sgn(jc,> — xj r ). 

Proof: The adjacency matrix of an undirected graph 
with no self-loops has (2) independent entries (for example, 
the upper triangular entries). Each one of these entries are 
a function of the vector of positions x. Hence, applying the 
chain rule, we have the following expansion for the partial 
derivative of the cost function with respect to the entries X{ r : 

df(x) . df da 



'pq 



l<p<q<n ^ a pq " X >> 



(7) 



Furthermore, a particular entry a pq depends solely on the 
position x p and x q ; therefore, we have that = if both 
p and q are different than i. Thus, for a fixed i, only the 
following summands in (|7]i survive 



dm 

dx ir 



y df da pi 
~l da pi dx ir 
^ df dajj 

dx ir 



y df da iq 
Pi da k dx ir 



pi da ij 
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where we have used that an = aji and a,, = 0. 

We now analyze each one of the partial derivatives in ([H} 
First, from |5]l and Q, we have that 
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where we have used Lemma [2] in the last equality. Second, 
from ([3} we have that 
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is not differen- 



for Xj r ^ Xj r (the fact that |ji 
tiable at x„ = Xj r does not affect our analysis). Let us 
define the antisymmetric matrix S r (x) = 



f ( r ) , 

y 



with 



s|' (1) =&ga(xi r — Xj r ). Hence, we can write 
Hadamard product as follows 

da: i(x) r , . , , , 

-MU. = -c\A(x)oS r (x)] {j . 

OXir 

Substituting Q and ( fTT) in ((8), we have 
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where both matrices A and S r depend on x. Then, from Q, 
we obtain the statement of our lemma. ■ 
An efficient relaxation of the spectral control Problem [3] 
results from controlling a truncated sequence of spectral 
moments (mi,...,m s ), for s < n. In this case, we can define 
a cost function 



* 1 2 

Mx) = £ ~{m k (x)-ml) , 
k=i * K 



(12) 



and an associated control law u — —V x f s (x). An explicit 
expression for u can be obtained by following the steps in 
the proof of Lemma [4] which result in 



w c 

U ir = 

n 



k=\ 
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(13) 



Although controlling a truncated sequence of moments 
is not mathematically equivalent to controlling the whole 
eigenvalue spectrum of the adjacency, we observe a very 
good overall matching between the eigenvalues obtained 
from the relaxed problem and the desired eigenvalue spec- 
trum, especially for the eigenvalues of largest magnitude, 
which are usually the most relevant in dynamical problems 
(Section [TvT>. 



B. Convergence Analysis 

To simplify convergence analysis of the closed-loop sys- 
tem pj, we restrict its dynamics to the open set & = 
{x I m k (x) > m* k , V 1 < k < s}^ We can ensure that x G & 
for all time by adding the barrier potential 



&,(*) = ! 



k=\ 



4k 



(m k (x) - 1 



(14) 



to |5]), for sufficiently small constants £1, . . . , E s > (assum- 
ing that the initial state of the system is already in J?). The 
convergence properties of the resulting closed loop system 
x= —V x (f(x)+b s (x)) are discussed in the following result. 

Theorem 5: Let {m£}| =1 denote a desired set of adjacency 
spectral moments and assume that x(Q) € & . Then, for 
sufficiently small £j . . . , £j > 0, the closed loop system 



(15) 



ensures that the moments {m/ ( (x)y k=l , where x = lim ; ^ TC x(f), 
approximate arbitrarily well the desired set {m?}i v 

2 It is shown in Theorem |5] that this construction does not restrict 
convergence of the state variables to the desired equilibria. 



Proof: The time derivative of f s (x) + b s (x) is given by 

b s {x)) = V x (f s (x)+b s (x))x 

= -||V«(/,(jt)+ft,(*))||i<0(16) 



d Jt im 



and so the closed loop system is stable and x will converge 
to a minimum of f s (x) +b s (x). Since b s (x) —> °° whenever 
x -> where = {x\m k (x) = m* k , V 1 < k < s} — 
{x | f s (x) = 0} denotes the boundary of the set & , equation 
( fT6| i also implies that the set & is an invariant of motion for 
the system under consideration. Let 

5*0 = {x | Xir(0) < Xj r (0) *,v < xj r , V /,;>} 

denote the polytope defined by the relative positions of the 
robots with respect to their initial configuration x(0), so that 
x(Q) € l~l & (see Fig. 1 I. In what follows we show that 
as £[,..., e s —> 0, the state variable x asymptotically reaches 
a value in the set £?o<~)dJF, where f s (x) = 0. 
To see this, observe first that 

lim (f s (x)+b s (x))=f s (x) 

£l,...,£,->0 

for all x € H i.e., for small enough E\ , . . . , £ s > 0, the 
potential / s (x) +b s (x) can be approximated by f s (x) in 
j£\ Moreover, note that f s (x) is convex in the set ^oHjF. 
Convexity of f s (x) follows from the fact that — c\xi r — Xj r \ is 
affine in f° r an y c > and, therefore, — Lc|x,y — Xi r \ is 
also affine for any number of terms in the summation. This 
implies that e~^ c \ x ' r ~ x i r \ is convex in as a composition 
of a convex and an affine function and, therefore, 



1 



trA k (x)-m* 



in 



-C\X ir —X\, 



■m k 



m k 



is also convex as a sum of convex functions. Clearly, 
j\rA k (x) — mi is nonnegative for any x <E ^oH,^. Since any 
power greater than one of a nonnegative and convex function 
is also convex, every one of the terms (~trA k (x) —fnfy , for 
k = 1, . . . ,s, is convex, which implies that f s (x) is also convex 
in l~l & . Taking the limit of ( fT6) > we have that 

d 



lim 



■ ,e, dt 



(fs(x)+b s (x))=f s (x)<0 



in ^qPiJ^. Therefore, convexity of f s (x) along with the 
condition f s (x) < implies that x will converge to a global 
minimum of f s (x) in ^qD^ (recall that & is an invariant 
of motion for the system under consideration). Since 

we conclude that the system will converge to a network with 
spectral moments {m k (x)}t =l that are almost equal to the 
desired {trit}l_ v The quality of the approximation depends 
on how small the constants £i , . . . , e s > are. 

What remains is to show that d& C For this, assume 
that d& <f_ Then, if the desired set of moments 



Fig. 1. Plot of the function e >*> r "i'K Note the boundary Xj r = xj r of two 
adjacent convex polytopes 8P, 



is realizable, there exists another polytope 3?* ^ such 
that d,j£ C 3 d *. Equivalently, there exists a configuration 
x* € such that A(x*) has eigenvalues {A,(x*)}" =1 with 
E" =1 Af(;t*) = m\ for all k = l,...,s. The result follows 
from the observation that ZPq can be obtained from 8?* 
by changing the relative positions of the robots in W 1 . 
Mathematically, this means that there exists a permutation 
matrix II e ^ndx.nd^ j e ; an orthogonal matrix with or 1 
entries, such that that if 3** = {x \ Ex > 0} with E £ W ldxnd , 
then 3* {) = {x\ E(Ux) > 0}. Therefore, if x* e 3 s *, then 
H T x* € ^o- Moreover, 

A(nV) ®/ d = n(A(x*) <g>/ d )n r , 

i.e., a permutation of the robots' coordinates results in a 
permutation of the entries of (A (x) where ® indicates 

the Kronecker product between matrices. Since n is orthog- 
onal, A(Tl T x*) and A(x*) ®I c i have the same eigenvalues. 
Therefore, there exists a configuration x = Yl T x* 6 sucn 
that A (I) has eigenvalues {A,-(jc)}" =1 with £" =1 ^f{x) =nit 
for all k = 1 , . . . , This implies that <3 J? C 3*$. 

In the above discussion, we have used the fact that if 
{y,}" =1 are the eigenvalues of Y 6 M" x " and {zi}'" =l are 
the eigenvalues of Z e R mxm , then the eigenvalues of Y ®Z 
are {y\z\ ,yiz m , ■ ■ ■ >y n Z\ ,y„z m }. This implies that the 
eigenvalues of A(x)(g)/</ are essentially the eigenvalues of 
A(x), each one with multiplicity d. ■ 

Remark 6 (Barrier functions): The barrier functions b s (x) 
are necessary in the proof of Theorem [5] If not there, the 
quantities m k (x) — m* are not guaranteed to be positive and, 
therefore, the potential f s (x) is not necessarily convex in 
3*0. This causes technical difficulties in ensuring a global 
minimum of f s (x). Also, the analysis in Theorem [5] assumes 
that e\,...,E s 0. This essentially restricts the influence of 
the barrier potential b s (x) to a small neighborhood of the 
set d^, so that the potential f s (x) remains unaffected out- 



3 The polytope &o essentially defines an ordering of the state variables. 



Equivalently, this means that the two graphs are isomorphic. 




Fig. 2. Hexagonal formation of mobile robots. 



side this neighborhood. In practice, the closer the constants 
El,..., e s are to zero, the better the approximation of the 
desired spectral moments. This approximation can be made 
arbitrarily good. 

Remark 7 (Distance metric): In the preceding analysis, 
we have employed the l\ norm (z = 1) as a distance metric 
to define the entries of the adjacency matrix (see ([3]l). This 
choice is mainly due to technical reasons, since it ensures 
that txA k (x) and, therefore, f s (x), are convex functions. From 
a practical point of view, the l\ metric has potential appli- 
cations in indoor navigation where the presence of obstacles 
forces the signals to propagate in grid-like environments. 
Nevertheless, our numerical simulations indicate that the 
control law herein proposed also performs well for z = 2 
(£2 norm). We leave a rigorous proof of this case for future 
work. 

IV. Numerical Simulations 

In this section we present examples of spectral design 
of mobile robot networks and discuss performance of our 
proposed approach. 

Example 8 (Hexagonal Formation): Consider the prob- 
lem of controlling the structure of a network of mobile robots 
to match the eigenvalue spectrum of the hexagonal network 
on 7 nodes shown in Fig. [2] The target eigenvalues of the 
weighted adjacency matrix of this formation are 

{l*} 7 i=1 ={-0.51, -0.47, -0.40, -0.40, 0.05, 0.05, 1.70}, 

and the corresponding spectral moments are 

{m* k } 7 k=l = {0, 0.53, 0.64, 1.22, 2.02, 3.47, 5.90}. 

The initial configuration of the mobile robot network is 
that of a random geometric graph on 7 nodes uniformly 
distributed in the square [0,1] x [0,1]. We applied the pro- 
posed control law ( fT5] l and studied the evolution of the 
spectral moments of the network's adjacency matrix. Fig. [3] 
shows the evolution of the second, third and fourth momentjj 

5 A similar behavior is observed for the higher order moments. Note that 
the first spectral moment is always equal to zero for simple graphs. 
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Fig. 3. Evolution of second (blue), third (green) and fourth (red) spectral 
moments towards the moments of the hexagonal formation. 



As expected, they all converge to the desired values. The 
asymptotic values of all spectral moments of the network 
are 

{m k }l =l = {0, 0.53, 0.65, 1.23, 2.04, 3.50, 5.95}, 

which are very close to the desired sequence of moments 
and slightly on the larger side, as predicted by 
Theorem [5] The eigenvalues of the weighted adjacency 
matrix of the final configuration are 

{A,}J =1 ={-0.52, -0.48, -0.42, -0.40, 0.02, 0.10, 1.70}, 

which are also very close to the desired values. Notice that 
our approach fits better those eigenvalues further away from 
the origin, since they are more heavily weighted in the 
expression of spectral moments, than those close to zero. 
An alternative approach to overcome this limitation would 
be to modify our cost function |5]l to assign more weight to 
eigenvalues of small magnitude. 

As discussed in Section |III-A| we also consider a relax- 
ation of Problem [3] that involves a truncated sequence of 
the first four moments of the network. In this case, the 



cost function is given by (12i and the closed loop system 



is defined in ( fT5j ), for 5 = 4. The asymptotic values of the 
first four spectral moments of the network are 

{m k } 4 k=l = {0, 0.55, 0.65, 1.26} 

and eigenvalues of the final configuration are 

{A. ; }J =1 = { 0.55, -0.50, -0.42, -0.31, -0.09, 0.18, 1.71}. 

which are remarkably close to the set of desired eigenvalues, 
especially those far from zero. The main advantages of 
using a relaxation of Problem [3] are (J) that it reduces the 
computational cost of the controller, since it only requires 
the terms trA k and A* -1 , and (ii) that it improves the stability 
of the numerical behavior of the gradient descent algorithm. 

Example 9 (Random Geometric Network): In this exam- 
ple we consider the problem of matching the eigenvalue 
spectrum of a particular realization of a random geometric 




graph on 10 nodes that are uniformly distributed in [0, 1] x 
[0,1]. This target realization is illustrated in Fig. |4j where 
the thickness of each edge is proportional to its weight. The 
eigenvalues of the weighted adjacency matrix of this network 
are 

{V>i=i = {-0.89, -0.85, -0.84, -0.79, -0.77, -0.68, 
-0.61, 0.02, 0.27, 5.16}, 

and the corresponding sequence of spectral moments is 

K} i>0 = {0, 3.11, 13.45,71.60, 368.36, 1905,...}. 

We applied the control law ( fT5j ) using the cost function 
(|5]l and starting with a different realization of the random 
geometric graph on 10 nodes in the square [0,1] x [0,1]. 
In this case, matching the whole set of spectral moments 
directly provides us with a very good approximation of the 
eigenvalues spectrum. As before, we also considered the 
relaxation of Problem [3] involving only a truncated sequence 
consisting of the first 4 moments. The evolution of the 
second, third and fourth moments of the network is shown 
in Fig. [5] The asymptotic values of the spectral moments of 
the network are 

{m k } k>0 = {0,3.13, 13.46,71.63,...}, 

which are remarkably close to the desired ones and, as in 
the previous example, slightly on the larger side. Similarly, 
the final set of eigenvalues is 

W}]=i = {-0-94, -0.87, -0.85, -0.83, -0.74, -0.59, 
-0.55, -0.27, 0.50, 5.16}, 

which is a very good fit to the given eigenvalues, especially 
for those far from zero. 

V. Conclusions and Future Research 

In this paper, we proposed a novel control framework to 
modify the structure of a mobile robot network in order to 
control the eigenvalue spectrum of its adjacency matrix. We 
introduced a novel abstraction of the eigenvalue spectrum by 



means of the spectral moments and derived explicit gradient 
descent motion controllers for the robots to obtain a network 
with the desired set of moments. Since the eigenvalue 
spectrum is uniquely determined by the associated spectral 
moments, our approach provides a way of controlling the 
eigenvalues of mobile networks. Convergence to the desired 
moments was always guaranteed due to convexity of the 
proposed cost functions. Efficiency of our approach was 
illustrated in nontrivial computer simulations. The adjacency 
matrix eigenvalue spectrum is relevant to the performance of 
many distributed coordination algorithms run over a network. 
Therefore, our approach is particularly useful in providing 
network structures that are optimal with respect to networked 
coordination objectives. 

Future work involves theoretical guarantees for the Eu- 
clidean distance metric (z = 2 in ([3])), as well as extension 
of our results to the spectrum of Laplacian matrix of the 
network. Moreover, the relaxation of Problem [3] to a trun- 
cated sequence of moments does not guarantee (mathemati- 
cally) a good fit of the complete distribution of eigenvalues. 
Therefore, a natural question is to characterize the set of 
graphs most of whose spectral information is contained in a 
relatively small set of low-order moments. 
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